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Executive  Summary 


This  report  presents  the  sequential  test  method.  The 
methods  are  described , a plan  and  computer  program,  to 
speed  the  use  of  the  method,  are  included.  This  testing 
method  could,  if  implemented  on  applicable  testing  pro- 
blems, reduce  test  time  and  costs  of  testing.  The  added 
benefit  of  reduced  energy  needs  are  inherent  in  this 
testing  method. 

The  text  was  originally  released  by  the  authors  in  1972. 
The  text  has  been  thoroughly  reviewed  and  is  applicable 
to  many  Tank-Automotive  system  and  component  testing  re- 
quirements. The  principal  efforts  involved  in  preparing 
this  report  were  performed  by  Mr.  John  Schmuhl,  who  is 
currently  employed  by  TARADCOM. 


1 • Introduction  and  Summary 

The  purpose  of  this  report  is  to  present  a method  and  a 
plan  which  could,  if  implemented,  reduce  the  amount  of  time 
needed  to  life  test  items  and  at  the  same  time  reduce  the 
costs  of  testing.  The  advantage  of  this  method  is  that  it  would 
not  be  necessary  to  test  an  item  through  the  complete  length  of 
time  specified  for  a test.  Necessary  decision  critiera  as  to  the 
acceptability  or  unacceptability  of  a test  can  be  determined 
much  earlier  thus  saving  time  and  money.  The  method  employed  to 
do  this  is  known  as  sequential  analysis  and  has  been  developed 
and  used  for  approximately  25  years.  A bibliography  of  articles 
and  related  topics  to  sequential  analysis  and  testing  is  given 
in  Appendix  IV. 


II . Sequential  Analysis 

One  hundred  percent  conclusive  and  valid  reliability  demon- 
strations usually  are  extremely  expensive  in  terms  of  time  and 
money.  To  prove  endurance  capabilities,  either  long  life  tests 
must  be  run,  or  if  time  is  a critical  factor,  large  numbers  o'f 
probably  expensive  items  must  be  placed  on  test  simultaneously. 

Sequential  analysis  takes  advantage  of  test  information  as 
it  is  accumulated  and  allows  for  previously  agreed  upon  decisions 
to  be  made  as  the  test  develops.  As  the  test  progresses  each 
failure  is  reported  and  plotted  on  a sequential  analysis  chart, 
such  as  that  shown  in  Figure  1.  This  plot  indicates  one  of  three 
possible  decisions  which  can  be  made  each  time  a failure  occurs: 
(1)  Reject  the  item,  either  individually  or  the  entire  lot  from 
which  it  was  drawn;  (2)  accept  the  item;  or  (3)  continue  the  test 
until  more  data  becomes  available. 

To  reach  these  decisions,  four  values  are  established  as 
criteria  against  which  test  results  are  compared  to  establish 
compliance  to  the  reliability  requirements  of  the  item.  These 
are  the  lowest  acceptable  mean  time  between  failure  (MTBF) 
usually  designated  "0X" , the  desired,  or  upper  limit  of  the  MTBF 
usually  designated  "60".  Associated  with  the  acceptance  of 
either  one  of  these  MTBF  values  are  risks,  normally  known  as 
"consumer's  risk"  and  "producer's  risk."  In' this  context 
consumers  risk  (statistically  known  as  Type  II  error  and  designa- 
ted "3")  is  the  probability  of  accepting  an  item  if  its  MTBF  is 
equal  to  9],  or  if  the  lot  is  actually  bad.  Producer's  risk 
(statistically  known  as  Type  I error  and  designated  "a")  is  the 
probability  of  rejecting  an  item  if  its  MTBF  is  equal  to  0O , or 
if  the  lot  is  actually  good.  This  assumes  that  0O  is  greater 
than  0i. 
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Four  values  are  required  in  light  of  the  different  interests 
of  the  "producer"  and  "consumer".  Because  the  true  MTBF  of  an 
item  cannot  be  determined  by  any  reasonable  sample,  two  different 
MTBF  values  are  determined,  one  high  and  the  other  low,  between 
which  the  "true"  value  is  expected  to  fall.  These  upper  (reject) 
and  lower  (accept)  limits  can  be  expressed  mathematically  and 
factored  into  the  sequential  analysis  chart.  The  most  common 
distribution  associated  with  sequential  life  testing  is  the 
exponential  distribution,  although  it  has  been  applied  to  the 
binomial,  Weibull  and  Poisson  distributions.  The  discussion  which 
follows  will  deal  with  the  exponential  distribution  as  most  work 
is  centered  on  it  and  it  is  mathematically  the  most  easy  to  deal 
with. 


III.  Mathematics  of  Sequential  Testing 

The  analysis  which  follows  assumes  the  following  situation: 
n items  are  placed  on  life  test  and  allowed  to  run  until  a failure 
occurs.  That  item  which  fails  is  either  replaced  or  repaired  and 
put  back  into  service,  thus  at  the  end  of  any  specified  period  of 
time  there  will  be  exactly  n items  on  test. . This  case  is  known  as 
the  Replacement  Case  as  opposed  to  the  Non-Replacement  Case,  wherein 
the  failed  items  are  not  replaced  or  repaired.  Mathematically, 
the  replacement  case  is  easiest  to  use  and  in  many  applications 
the  most  reasonable.  In  addition,  the  underlying  life  distribu- 
tion will  be  assumed  to  be  exponential  with  probability  density 
function: 

1 ”x/0 

f(X,6)  = ie  ; X > 0.  . 

The  variable  "x"  represents  time  and  the  unknown  parameter 
0(>O)  can  be  thought  of  physically  as  the  mean  life.  What  is 
done  in  sequential  analysis  is  to  test  the  simple  hypothesis 
that  the  true  MTBF,  9,  equals  9q  (the  upper  limit  of  the  MTBF) 
against  the  simple  alternative  hypothesis  that  9 equals  9i  (the 
lower  limit  of  the  MTBF) . The  test  is  carried  out  by  drawing 
n items  at  random  from  the  population  and  placing  them  all  on 
life  test.  The  basic  rationale  for  the  formulation  of  sequential 
tests  is  derived  from  work  done  by  A.  Wald  [44]..  Wald’s  work  on 
sequential  analysis  can  be  used  virtually  without  modification 
in  a situation  where  decisions  are  made  continuously.  A word  on 
notation  is  in  order  at  this  point.  When  referring  to  the 
"hypothesis  that  the  true  MTBF  equals  &£ , it  is  common  to 
abbreviate  this  as: 

Ht  : 6 « 0j 

and  the  "hypothesis  that  the  true  MTBF  equals  6o"  as: 

Ho  : 9 = 0Q 
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We  v/ish  to  test  Hg  s 9 * 0g  against  Hj  : Q » 0j  with  producer's 
and  consumer's  risk  a and  6,  respectively.  Since  information  is 
available  continuously,  a continuous  analogue  of  the  sequential 
probability  ratio  test  of  Wald  can  be  used.  The  decision  to 
accept,  reject,  or  continue  the  test  depends  on: 

3 < (6o/0i)r  exp  [-(1/0!  - 1/0O ) V (t) ] < A (1) 

where  3 and  A are  constants,  depending  on  a and  & such  that 
B < 1 < A and 


B = 


3 

Ta^aT 


a = 4CS2 
a 


(la, lb) 


The  variable  "r"  represents  the  number  of  failures  up  to  time  "t". 
The  decision  to  continue  testing  is  made  as  long  as  the  inequality 
(1)  holds.  At  the  time  the  experiment  is  stopped,  if  the  first 
inequality  in  (1)  is  violated  Hg  is  accepted;  if  the  second 
inequality  is  violated  Hi  is  accepted.  V(t)  is  a statistic  which 
can  be  interpreted  as  the  total  life  observed  up  to  time  t.  In 
the  replacement  case: 

V(t)  = nt  (2) 

Inequality  (1)  can  be  reformulated  more  conveniently  as: 

-hi  + rs  < V (t)  < hg  + rs  (3) 


where  hj , hg,  and  s are  positive  constants  given  by: 


^0 


hx 


s 


r 


-In  (B) 

1/0^  1/eg 

In  (A) 

i/ei~i/e0 

In  ( 6 o/9 l ) 

1/  9 g 

failure  number 


In  plotting  the  sequential  test  chart,  the  first  part  of  inequality 
(3)  forms  a straight  line  representing  the  reject  boundary  while 
the  second  part  of  inequality  (3)  forms  a straight  line  represent- 
ing the  accept  boundary.  Graphically,  this  is  illustrated  in 
Figure  2. 
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SEQUENTIAL  TEST  CRITERIA 


Further,  it  can  be  shown  that  the  probability  of  accepting 
Ho,  L ( 8 ) , when  8 is  the  true  parameter  value,  is  given  approximately 
by  a pair  of  parametric  equations: 


L(0) 


(7a,  7b) 


e = (eo/ei)h-i 
h(i/ex-i/e0) 

by  letting  the  parameter  h run  through  all  real  values.  The  values 
of  L(8)  at  five  points  8 = 0,  0i,  S,  0^,  and  00  are  enough  to  sketch 
the  entire  curve  (known  as  the  "Operating  Characteristics"  curve) . 
These  values  are,  respectively,  0,  0,  In (A) /In (A) -In (B) , 1-a,  and  1. 


Additionally,  it  is  possible  to  determine  the  expected  number 
of  observations  required  to  reach  a decision  when  6 is  the  true 
parameter.  This  quantity  is  abbreviated  E0  (r)  and  is  given 
approximately  by: 

( L (8 ) InS  + [ 1-L  (0)1  InA  _ ^-1,(8)  (hQ+hj) 

TnToJTeJT  - e (1/ei-i/0o)  s-e  ?e  ^ s <9a) 


E0  \ 


-In (A) In (B) 
f in  ( © 0/ 0'jL ) ] 


= hQ£L 
sz 


;e  = s (9b) 


Letting  k = 0o/0n,  three  important  values  of  E0  (r)  become  parti- 
cularly simple  when  0=6  3.,  s,  or  0q.  They  are: 

E01  (r)  * [ 01n(B)  + (1-0)  In  (A)  ] / [In  (k)  - (k-1)  /k] 

Es(r)  - -In (A) In (B)/[ln (k) 3 2 


Ee0  (r)  = [ (1-a) In (B)  + aln (A) / [ln(k)  - (k-1)] 


(10a) 

(10b) 

(10c) 


In  the  replacement  case  where  the  number  of  items  on  test 
throughout  is  the  same,  namely  n,  it  can  be  shown  that  the  expected 
waiting  time,  Eq  (t) , before  a decision  is  reached  if  0 is  the  true 
MTBF , is  given  by: 

Ee  (t)  = (8/n) E0  (r ) 


It  is  prudent,  at  this  point,  to  stop  and  summarize  the 
sequential  decision  criteria: 


If  -hx+ rs  <V(t)<  h0+  rs,  continue  the  test,  (12a) 
If  V (t)  > hQ  + rs,  stop  the  test  and  accept  Hq,  ■ (12b) 
If  V (t)  _<  -hj  + rs,  stop  the  test  and  reject  Hq,  (12c) 


7 


It  can  be  seen  that  it  is  possible  to  satisfy  (12a)  for  all  t, 
hence  making  it  impossible  to  ever  make  an  accept  or  reject 
decision.  In  other  words,  the  life  test  could  run  forever.  What 
is  done  to  prevent  this  situation  from  occurring  is  either  to  set 
a maximum  time,  Tq,  at  which  to  terminate  the  test  or  a maximum 
number  of  failures,  ro,  at  which  to  terminate  the  test  or  both. 

This  guarantees  that  the  test  will  run  no  longer  than  time  To  or 
the  time  at  which  rg  failures  occur.  Usually  one  knows  Tq  and  must 
find  Tq.  This  is  easily  done,  provided  one  has  a table  of  chi- 
square  values  (Appendix  II) . Epstein  [17]  has  shown  that  ro  is  the 
smallest  integer  for  which: 


X l-g,2r 


S,2r 


This  can  be  easily  related  to  Tq  from  Epstein  [17]  by: 


(13) 


e0X2i-a,2r 

T = 1 

0 2n 


(14) 


solving  this  for  6q  yields: 


2n  Tn 

Y2 

* l-g,2r 


and  substituting  this  into  (13)  gives: 


l-g,2r 


61X  l-g , 2r 
2nT, 


X B/2r  - —0 

Rearrange  (16)  and  simplify: 


iSlfl  > 
01  - 


e,2r 


(15) 


(16) 


(17) 


ro  = r now  the  largest  integer  such  that  (17)  is  true.  To  and  ro 
are  commonly  referred  to  as  the  truncation  time  and  truncation 
failures,  respectively. 

This  essentially,  summarizes  the  basics  of  exponential 
sequential  analysis  as  applied  to  life  testing.  For  more 
detailed  presentations  of  this  subject  the  bibliography  in  Appendik 
IV  should  be  consulted. 
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Examples  of  Sequential  Analysis  in  Lifa  Testing 

EXAMPLE  1*  Consider  a sample  of  one  item  which  is  placed 
on  replacement  life  test.  The  underlying  distribution  is  expon- 
ential with  the  true  mean  time  between  failure  located  somewhere 
between  320.0  hours  and  500.0  hours.  A consumer's  and  producer's 
risk  of  20%  is  decided  upon  and  a truncated  sequential  test  plan 
is  desired. 

SOLUTION: 


0O  = 500.0  hours  ot  = 0.20 

0j  = 320.0  hours  B = 0.20 

n = 1 


To  determine  rg , the  truncation  failures,  (13)  should  be  used 
along  with  Appendix  (II) . From  Appendix  II  it  is  seen  that  when 
v = 2r  = 30: 

X2 i -a , 2r  xZ 0,80, 30 

x2  ~ x2 

0,2r  0,20,30 

Thus : 

rQ  = 2r/2 

= 15 


23.364 
= 36.250 


= 0.644  > 


- 0, 


320.0  _ n 


Knowing  rn  it  is  possible  to  determine  the  truncation  time  from 


(14)  : 


®0X" 


To  = 


i-g,2r 


2n 


500.0x2q.80  , 30 
2(1) 


= 500.0  (23 . 364 ) /2 

= 5840.0  hours 


The  sequential  decision  criteria  are  developed  from  (4) , (5) , 
and  (6) , namely: 


B _ 0.20 

(l-o)  0.80 


(1-B)  = 0.80 

a 0.20 


4 


h0  = -In  B 

1/e1-i/e0 

= 1233.5 


-In  (1/4) 
1/320  - 1/500 
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In  A 


1 1/0  x “ 1/6 0 
= 1233.5 

s * in  (60/01) 

l/e i - i/e0 
= 396.6 


_ In  (4) 

1/320  - 1/500 

= In  (500/320) 
1/320  - 1/500 


Thus: 

If  -1233.5  + 396. 6r  <t<  1233.5  + 396.6r,  continue  the  test. 


If  t > 1233.5  + 396. 6r,  stop  the  test  and  accept  H0, 
If  t £ 1233.5  + 396. 6r,  stop  the  test  and  reject  Hc. 


Graphically  these  decision  rules  appear  in  Figure  3.  Assume  the 
following  failure  times  were  noted  during  the  test  (in  hours) : 

400,  1000,  1200,  2400,  3400 

Note : These  are  measured  from  time  equal  to  zero;  they  are  not 
time  between  failure. 

These  points  are  plotted  on  Figure  3 and  shown  in  Figure  4. 
It  can  be  seen  that  the  test  can  be  terminated  after  2800  hours, 
before  the  fifth  failure  occurs,  with  a decision  to  accept  the 
hypothesis  that  the  true  mean  time  between  failure  is  500  hours 
(or  greater) . It  is  interesting  to  note  that  had  no  failure 
occurred  before  the  first  1233.5  hours,  the  test  could  have  been 
terminated  with  an  accept  decision. 

Points  for  the  Operating  Characteristics  curve,  expected 
failure,  and  expected  time  curves  are  summarized  below: 


0 

L (0 ) 

E0  (r) 

E0  (t) 

0.0 

0.0 

3.1 

0.0 

320.0 

0.2 

9.7 

3092.2 

396.6. 

0.5 

9.7 

3836.8 

499.7 

0.8 

7.2 

3586.8 

00 

1.0 

0.0 

1233.5 

These  are  plotted  in  Figures  5,  6,  and  7 respectively. 

In  order  to  simplify  the  calculations  required,  a computer 
program  has  been  written  and  implemented  to  develop  sequential 
test  plans.  The  program  written  in  FORTRAN  IV  language  has  been 
developed  for  the  GE-440  time-sharing  system.  The  program  is 
interactive  and  a listing  of  it  is  given  in  Appendix  III.  In 
order  to  illustrate  its  use  the  following  example  is  given. 
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TERMINATE  TEST  AT  15  FAILURES  WITH 
REJECT  DECISION,  TERMINATE  TEST  5840 
HOURS  WITH  ACCEPT  DECISION 


(j)  saamivj 
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( e)  t: 
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EXAMPLE  2:  Consider  an  engine  which  is  to  be  tested  for 
a maximum  of  2000.0  hours.  In  order  to  accept  the  engine  it 
must  demonstrate  a minimum  mean  time  between  failure  of  151.0  hours. 
A producer's  risk  of  10  percent  is  desired  with  only  a 1 percent 
consumer's  risk.  A sequential  testing  plan  is  required  assuming  the 
exponential  distribution. 

SOLUTION: 


The  following  data  is  given: 

T0  = 2000.0  hours  a ” 0.10 

= 151.0  hours  3 - 0.01 

n = 1 


The  program,  as  written,  will  compute  0O  and  rD,  however,  they 
will  be  manually  calculated  in  order  to  demonstrate  the  method. 
This  can  be  done  using  (17) : 

2nTo  = 2 (1) (2000.0)  _ oc  c 
6X  151.0 


rD  = r is  the  largest  integer  such  that: 


26.5  > x: 


3,2r 


= X‘ 


0.01, 2r 


When  2r  = 12,  Xq.01,  12  “ 26.217,  therefore, 


r0  = 6 


Having  determined  r^  it  is  possible  to  obtain  0o  from  (15): 


2nT 


fl. 


* l-a,2r 


2 (1) (2000.0) 

X2  0 . 90 , 12 


4000 

6.304 


= 634.1  hours 


The  sequential  decision  criteria  can  be  found  exactly  as  before.  4 
They  are  given  below: 

If  ~454.4  + 294. 4r  < t < 891.8  + 284. 4r,  continue  the  test. 


If  t > 891.8  + 284. 4r,  stop  the  test  and  accept  Hgf 

If  t _<  -454.4  + 284. 4r,  stop  the  test  and  reject  Hq. 

Graphically,  these  decision  criteria  are  shown  in  Figure  8. 

It  is  seen  that  the  minimum  time  that  the  test  could  run  is 
891.8  hours  for  an  accept  decision  or  the  time  at  which  the 
sixth  failure  occurs  for  a reject  decision.  The  complete  computer 
dialogue  and  output  for  this  example  is  shown  in  Figure  9 . In 
addition,  the  program  can  automatically  vary  a and  6 between 
four  values,  namely,  0.01,  0.05,  0.10,  0.20  so  that  16  possible 
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TERMINATE  TEST  AT  6 FAILURES  WITH 
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FIGURE  8 - SEQUENTIAL  TEST  CHART  FOR  EXAMPLE 


RUM 


SEQTST  15:25  ACTS  06/27/72  TUE. 


ENTER  N ,THET A 1 ,T0 
INPUT:? 0146 
? 1,151. ,2000. 


ENTER  0 IF  FULL  DISPLAY  IS  DESIRED 
ENTER  1 IF  NOT  DESIRED 
INPUT : BE  1 54 

? 1 


ENTER  ALPHA, BETA  DESIRED 
INPUT :00163 
? 0.10,3.01 


SEQUENTIAL  TESTING  PLANS 
EXPONENTIAL  DISTRIBUTION 
(WITH  REPLACEMENT) 
**#**#***##*#■***#*#***#* 


ALPHA  - 0.10 
BETA  « 0.01 
N » 1 
T0  - 2000.0 
THETA1  - 151.0 

****#«**##*###*# 


R ■ 6 

THETA0  - 634.1 

a######*###**#** 

H0  - 891. G 

HI  - 454.4 

S - 284.4 


FIGURE  9 


1$ 


plotting  points 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

891.  B 

1 

0.0 

1176.2 

2 

114.4 

1460.6 

3 

398.8 

1745.0 

4 

6C3.2 

2000,0 

5 

967.6 

2000.0 

6 

1252.0 

2000.0 

************************************************************** 


THETA 

L (THETA ) 

E (THETA  ,R  ) 

E (THETA  ,T ) 

0.0 

0.000 

1.6 

0.0 

151.0 

0.010 

3.3 

499.1 

284.4 

0.338 

5.0 

1424.9 

634.1 

0.900 

2.2 

1373.0 

INFINITY 

1.000 

0.0 

891.8 

********** **************************************************** 

STOP 

RUNNING  TIME:  4.4  CPUS  ELAPSED  TIME : 5.0  CPUS 


FIGURE  9 (Continued) 
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sequential  test  plans  can  be  generated.  The  output  for  these  are 
given  in  Figures  10  through  25. 

V.  Conclusions 


This  report  has  dealt  with  the  subject  of  sequential  analysis 
as  applied  to  life  testing.  Sequential  testing  has  as  its  main 
advantage  the  capability  of  reducing  the  amount  of  test  time 
required  in  order  to  make  a decision  regarding  the  acceptability 
or  unacceptability  of  an  item  or  its  components.  This  can  mean 
reduced  costs  and  the  possibility  for  a larger  number  of  tests. 

From  a statistical  standpoint  sequential  testing  is  sound  and  in 
practice  could  be  administered  simply.  The  remainder  of  this  report 
is  devoted  to  Appendices  and  a bibliography  which  can  be  consulted  if 
further  investigation  into  the  subject  is  necessary.  It  is  this 
author's  opinion  that  references  [17]  and  [20],  together,  provide 
the  most  detailed  and  complete  treatment  of  sequential  testing 
and  was  used  considerably  in  the  preparation  of  this  report. 


SEQUENTIAL  TESTING  PLANS 

EXPONENTIAL  DISTRIBUTION 

(WITH  REPLACEMENT) 

#*######*######*#*##**## 

ALPHA  - 0.01 

BETA  - 0.01 

N - 1 

T0  - 2000.0 

THETA  1 - 151.0 

**#*#«***#**##** 

R ■ 6 

THETA0  - 1134.3 

**************** 

H0  m 800.4 

HI  - 8 00.*  4 

S - 351.2 

PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

800  .-4 

1 

0.0 

1151.7 

2 

0.0 

1502.9 

3 

253.3 

1B54.2 

4 

604.6 

2000.0 

5 

955. B 

2000.0 

6 

1307. 1 

2000.0 

* **********#■**#*****##**#*#*##**#******#*****#*  #****#•**#**#**# 


THETA 

L(THETA) 

E (THETA  ,R ) 

E(THETA.T) 

0.0 

0.000 

2.3 

0.0 

151.0 

0.010 

3.9 

591.5 

351.2 

0.500 

5.2 

1824.0 

1134.3 

0.990 

1.0 

1 136.3 

INFINITY 

1.000 

0.0 

800.4 

•a*#*##*#*##****##*****#***#*##*#***#*#*#****##**##*******#*#*# 


FIGURE  10 


ALPHA  - 0.01 

BETA  - 0.05 

N - 1 

T0  - 2000.0 

THETA  1 - 151.0 

ft***###**##****# 

R - 8 

THETA0  - 692.9 

#**#*#**#*##**** 


H0  - 576.5 
HI  « 879.2 
S - 294.2 

PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

576.5 

1 

0.0 

870.6 

2 

0.0 

1164.8 

3 

3.3 

1459.0 

4 

297.5 

1753.2 

5 

591.6 

2000.0 

6 

885.8 

2000.0 

7 

1 180.0 

2000.0 

8 

1474.2 

2000.0 

**********#*****#*#*#*##*###*#***###*###***##*■#*##*##***##***# 


theta 

L(THETA) 

E (THETA ,R ) 

E (THETA, T 

0.0 

0.000 

.3.0 

0.0 

151.0 

0.050 

5.6 

850.5 

294.2 

0.604 

5.9 

1722.9 

692.9 

0.990 

1.4 

976.4 

INFINITY 

1.000 

0.0 

576.5 

*#*#*#*#*#**  #****##**###****■*•****#**  *####*#  ft####***###*#****** 


FIGURE  11 


22 


ALPHA  - 0.01 

BETA  - 0.10 

N - 1 

T 0 - 2000.0 

THETA  1 - 151.0 

**************** 

R 9 

THETA0  - 573.3 

**************** 

H0  « 470.0 

HI  - 922.4 

S » 273.5 

PLOTTING  POINTS 

FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

470.0 

1 

0.0 

743.4 

2 

0.0 

1016.9 

3 

0.0 

1290.4 

4 

171.5 

1563.9 

5 

445.0 

1837.4 

6 

718.5 

2000.0 

7 

992.0 

2000.0 

8 

1265.5 

2000.0 

9 

1539.0 

2000.0 

************************************************************** 


THETA 

L(THETA) 

E (THETA  ,R ) 

E (THETA  ,T ) 

0.0 

0.000 

3.4 

0.0 

151.0 

0.100 

6.4 

965.5 

273.5 

0.662 

5.8 

1585. 1 

573.3 

0.990 

1.5 

872.0 

INFINITY 

1.000 

0.0 

470.0 

******  *****************************************  *************** 


FIGURE  12 
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ALPHA  - 0.01 

BETA  - 0.20 

N « 1 

T0  « 2000.0 

THETA  1 - 151.0 

ft*****##**#*#*## 

R - 10 

THETA0  - 486.4 

««**#****##**#** 

H0  - 350,2 

HI  - 959.6 

S - 256.1 

PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

350.2 

1 

0.0 

606,4 

2 

0.0 

862,5 

3 

0.0 

1118.7 

4 

65.0 

1374.8 

5 

321.1 

1631.0 

6 

577,3 

1887. 1 

7 

833.4 

2000 . 0 

8 

. 1089.6 

2000.0 

9 

1345.7 

2000.0 

10 

1601.9 

2000.0 

*#*#*■*#**#*•*■#****###*#**##***##**#**##■*####***#******#****#*## 


THETA 

L(THETA) 

E(THETA,R) 

E (THETA ,T) 

0.0 

0.000 

3.7 

0.0 

— 

151.0 

0.200 

6.6 

1001.9 

256.1 

0.733 

5.  1 

1312. 1 

* 

486.4 

0.990 

1.5 

712.3 

- 

INFINITY 

1.000 

.0.0 

350.2 

#***#**«*# *«*«*## #«#*###*#*#*#*###**#*#*###*##*#*****«#* #«###* 


FIGURE  13 
2k 


ALPHA 

- 0.05 

BETA 

- 0.01 

N 

■ 1 

T 0 

- 2000.0 

THETA 

1 - 151.0 

**************** 

R 

« 6 

THETA0  - 766.6 

**************** 

H0 

- 856.3 

HI 

- 561.4 

8 

- 305.5 

PLOTTING  POINTS 

FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

856.3 

1 

0.0 

1161.8 

2 

49.6 

1467.3 

3 

355.1 

1772.8 

4 

660.6 

2000.0 

5 

966.1 

2000.0 

6 

1271.6 

2000.0 

************************************************************** 


THETA 

l(theta) 

E(THETA,R) 

E (THETA  ,T 

0.0 

0.000' 

1.8 

0.0 

151.0 

0.010 

3.5 

534.8 

305.5 

0.396 

5.2 

1573.6 

766.6 

0.950 

1.7 

1305.8 

infinity 

1.000 

0.0 

856.3 

************************************************************** 


FIGURE  14 
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0.05 
0.05 
. 1 

2000,0 

151.0 

###***#•**##*#*** 

R - 8 

THETA0  - 502.8 

#*#*#***#*##***# 


H0  - 635.4 
HI  - 635.4 
S - 259.6 

PLOTTING  POINTS 


ALPHA  - 

BETA 

N 

T 0 

THETA  1 - 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

635.4 

1 

0,0 

895.0 

2 

0.0 

1154.7 

3 

143.4 

1414.3 

4 

403.0 

1673.9 

5 

662.6 

1933.5 

6 

922.2 

2000.0 

7 

. 1181.8 

2000.0 

8 

1441.4 

2000.0 

********###*#*******#**#•■***■*•#**#*#***###■**##*#***•***##*•***#*#* 


THETA 

L(THETA) 

E(THETA.R) 

E(THETAtT 

0.0 

0.000 

2,4 

0.0 

151.0 

0.050 

5.3 

795.1 

259.6 

0.500 

6.0 

1555.4 

502.0 

0.950 

2.4 

1182.4 

INFINITY 

1.000 

0.0 

635.4 

***#*****##*#**#****##*###****#*#****#********#**#*##***#*#*#.* 


Figure  15 
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ALPHA  - 0.05 
BETA  - 0.10 
N - 1 
T 0 « 2000.0 
THETA1  - 151.0 

**************** 


R - 9 

THETA0  - 426.2 

**************** 


H0  - 526.5 
HI  - 675.9 
S - 242.7 

PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

526.5 

1 

0.0 

769.1 

2 

0.0 

1011.8 

3 

52.1 

1254.4 

4 

294.7 

1497.1 

5 

537.4 

1739.7 

6 

700.0 

1982.4 

7 

1022.7 

2000.0 

8 

1265.4 

2000.0 

9 

1500.0 

2000.0 

************************************************************** 


THETA 

L(THETA) 

E(THETA.R) 

E( THETA, T 

0.0 

0.000 

2.8 

0.0 

151.0 

0.100 

6.  1 

915.4 

242.7 

0.562 

6.0 

1466.4 

426.2 

0.950 

2.5 

1082.8 

INFINITY 

1.000 

0.0 

526.5 

********* ***************************************************** 


FIGURE  16 
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ALPHA  - £5.05 

BETA  - 0.20 

N - 1 

T0  - 2000.0 

THETA1  - 151.0 

*###*#***«****#* 

R - 10 

THETA0  - 36B.8 

#*####*«#*#***«# 


H0  - 398.4 

HI  - 708.9 

S » 228.3 


PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

398.4 

1 

0.0 

626.7 

2 

0.0 

855.0 

3 

0.0  . 

1083.4 

4 

204.4 

1311.7 

5 

432.7 

1540.0 

6 

661.1 

1768.3 

7 

889.4 

1996.7 

8 

1117.7 

2000.0 

9 

1346.0 

2000 . 0 

10 

1574.4 

2000. 0 

■#***********•**#■*■*■****###*#**  **•**###*#*■*#****#****#*  ##*****##** 


THETA 

L (THETA ) 

E (THETA, R) 

E (THETA, T ) 

0.0 

0.000 

3. 1 

0.0 

151.0 

0.200 

6.3 

951.9 

228.3 

0.640 

5.4 

1236.9 

368.8 

0.950 

2.4 

900.5 

INFINITY 

1.000 

0.0 

398.4 

###***# *****##*#*#*#**#*#*#**********####****##*# #*#**#####*#* 


FIGURE  17 


28 


ALPHA  - 0.10 
BETA  - 0.01 
N - 1 
T 0 - 2000.0 
THETA  1 - 151.0 

**************** 


R - 6 

THETA0  - 634.1 

*#*«****#*«**««* 


H0  - 891.8 
HI  - 454,4 
S « 284.4 

PLOTTING  POINTS 


fail.  no. 

REJECT 

ACCEPT 

0 

0.0 

891.8 

1 

0.0 

1176.2 

2 

1 14.4 

1460.6 

3 

398.8 

1745.0 

4 

683.2 

2000. 0 

5 

967.6 

2000.0 

6 

1252.0 

200  0.0 

»*«#««*»««**«««##*««*«*****«#**«««*f***#*«*********«*«#*«##*#* 


THETA 

L (THETA) 

E (THETA, R) 

E (THETA, T ) 

0.0 

0.0  00 

1.6 

0.0 

151.0 

0.010 

3.3 

499.1 

204  .4 

0.338 

5.0 

1424.9 

634 . 1 

0.900 

2.2 

1373.0 

INFINITY 

1.000 

0.0 

891.  B 

#**#****#**#***##***#«*##*#**##**#**#**#*#**#**##*##********#* 


FIGURE  18 


29 


ALPHA  - B.10 

BETA  - 0.05 

N «=  1 

T0  « 2000.0 

THETA  1 - 151.0 

******■****##*#■»* 

R - 8 

THETA0  - 429.3 

«*##*#*******#«* 


H0  - 673.2 

HI  - 524.4 

S - 243.4 


PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0..0 

673.2 

1 

0.0 

916.6 

2 

0.0 

1160.0 

3 

205.8 

1403.4 

4 

449.2 

1646.8 

5 

692.6 

1890.2 

6 

936.0 

2000.0 

7 

1179.4 

2000.0 

8 

1422.7 

2000.0 

•***•#***#**##■***#*****#*  #**#**#*###*##**##****####*##**•***#*##  * 


THETA 

L (THETA) 

E (THETA  ,R ) 

E (THETA ,T 

0.0 

0.000 

2.2 

0.0 

151.R 

0.050 

5.0 

759.2 

243.4 

0.438 

6.0 

1450.4 

429.3 

0.900 

3.0 

1277.9 

INFINITY 

1.000 

0.0 

673.2 

•*****#■*##***  ^**  ********************  #»##***«#**#*  «#**#**#*«#**# 


FIGURE  19 
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ALPHA  - 0,10 
BETA  - 0,10 
N - 1 
T 0 - 2000.0 
THETA  1 - 151.0 

**************** 


R . 9 

THETA0  - 368.0 

**************** 


H0  - 562.7 
HI  - 562.7 
S - 228.1 

PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

562.7 

1 

0.0 

790.8 

2 

0.0 

1018.9 

3 

121.7 

1247.0 

4 

349.8 

1475.1 

5 

577,9 

1703.2 

6 

806.0 

1931.3 

7 

1034,1 

2000.0 

8 

1262.2 

2000.0 

9 

1490.3 

2000.0 

*****  ****  *************************  *********  ******************  * 


THETA 

L(THETA) 

E(THETA,R) 

E (THETA  ,T 

0.0 

0.000 

2.5 

0.0 

151.0 

0,100 

5.e 

881. 4 

220.1 

0.500 

6.1 

1387.8 

360.0 

0.930 

3.2 

1 1C4. 1 

INFINITY 

1.000 

0.0 

562.7 

************************************************************** 


FIGURE  20 
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ALPHA  - 0.10 

BETA  - 0.20 

N - 1 

TC  - 2000.0 

THETA  1 - 151.0 

#*#*##**##•#***** 

R - 10 

THETA0  - 321.4 

**************** 


H0  - 420.4 

HI  - 592.3 

S - 215.1 


PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

428.4 

1 

0.0 

643.6 

2 

0.0 

858.7 

3 

53. 1 

1073.0 

4 

268.2 

1289.0 

5 

483.4 

1504.  1 

6 

698.5 

1719.2 

7 

913.6 

1934.4 

8 

1128.8 

2000.0 

9 

1343.9 

' 2000.0 

10 

1559.0 

2000.0 

*#**#*#■»#*#****#*###**#*****■»**##***##*****#*#**•#******##**### 


THETA  L(THETA)  E(THETA.R) 


E (THETA, T ) 


0.0 

151.0 

215.1 
321.4 

INFINITY 


0.000 

0.200 

0.560 

0.900 

1.000 


2.8 

0.0 

6.1 

913.9 

5.5 

1179.6 

3.1 

987.3 

0.0 

428.4 

**####**##***####*###########**######****#***##**#*#**#**#**#* 


FIGURE  21 


32 


ALPHA  - 0.20 

BETA  « 0.01 

N - 1 

T0  - 2000.0 

THETA  1 - 151.0 

*#«*«#«**«#«*#*# 

R - 6 

THETA0  - 511.5 

ft*************** 


H0  - 938.8 

HI  - 342.7 

S - 261.4 


PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

938.8 

1 

0.0 

1200.2 

2 

180.1 

1461.6 

3 

441.5 

1723.0 

4 

702.9 

1984.4 

5 

964.3 

2000.0 

6 

1225.7 

2000.0 

*******#«■**#»**#*#***#*#»##*****#*********#«•**##**#**#*♦#***#* 


THETA 

L(THETA) 

E(THETA,R) 

E(THETA.T) 

0.0 

0.,  0 00 

1.3 

0.0 

151.0 

0.010 

'3.0 

451.1 

261.4 

0.267 

•4.7 

1230.7 

511.5 

0.800 

2.7 

1395.9 

INFINITY 

1.000 

0.0 

938.8 

*#*  #*■#■*■*#*•*■*#**#■*#**#**#-#*#■*##**#*#•**#*#•*************##*  ##*#•** 


FIGURE  22 
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ALPHA  - P.20 

BETA  - 0.P5 

N ■ .1 

T 0 - 2000.0 

THETA  1 - 151.0 

**#■#**##**###*** 


R ■ 8 

THETA0  - 358.3 

#*###****#*•*•**** 

H0  - 723.6 

HI  - 406.7 

S - 225.5 

PLOTTING  POINTS 


PAIL.  NO. 

REJECT 

ACCEPT 

fi 

0.0 

723.6 

1 

0.0 

949.  1 

2 

44.4 

1174.7 

3 

269.9 

1400.2 

4 

495.4 

1625.7 

5 

720.9 

1851.2 

6 

946.5 

2000.0 

7 

1172.0 

2000.0 

8 

1397.5 

2000.0 

***************************** a**#**#**##****#**#* ************* 


THETA 

L(THETA) 

E (THETA ,R ) 

E (THETA ,T 

0.0 

0.000 

1.6 

0.0 

151.0 

0.050 

.4.7 

709.5 

225.5 

0.360 

5.8 

1304.8 

358.3 

0.800 

3.7 

1342.6 

INFINITY 

1.000 

0.0 

723.6 

*#**•##*#  #****■##**#  ##*###•*****#  #*#*******#****#**###  *#**#**•**#* 


FIGURE  23 


ALPHA  - 0.20 

BETA  - 0.10 

N * 1 

T 0 - 2000.0 

THETA1  - 151.0 

*##*#*  **#*#■**  *** 


R - 9 

THETA0  - 310. B 

#**#******#**### 


H0  - 610.6 
HI  - 441.7 
S « 212.0 

PLOTTING  POINTS 


FAIL.  NO. 

REJECT 

ACCEPT 

0 

0.0 

610.6 

1 

0.0 

822.6 

2 

0.0 

1034.6 

3 

194.4 

1246.7 

4 

406.4 

1450.7 

5 

61B.4 

1670.7 

6 

830.4 

1882.7 

7 

1042.4 

2000.0 

B 

1254.5 

2000.0 

9 

1466.5 

2000.0 

* * *******#•***#*##***###***#  *******#■**#*#*#**#**##*#***■#*  *****  * 


THETA 

L ( THETA ) 

E (THETA  ,B  ) 

E (THETA ,T 

P.0 

0.000 

2.1 

0.0 

151.0 

0.  100 

• 5.5 

832.6 

212.0 

0.420 

6.0 

1272.0 

310.0 

0.800 

4.0 

1258.6 

INFINITY 

1.000 

0.0 

610.6 

******************###******#**#*##***#*#**##******##********## 


FIGURE  24 


ALPHA  - 0.20 
BETA  - 0.20 
N 1 
T0  - 2000.0 
THETA 1 - 151.0 

**•*«*•*•***«*#* 


R - 10 

THETA0  - 274.2 

**************** 


H0  - 466.0 

HI  - 466.0 

8 - 200.5 


PLOTTING  POINTS 


FAIL 

. 1 

NO. 

REJECT 

ACCEPT 

0 

0.0 

466.0 

1 

0.0 

666.4 

2 

0.0 

866.9 

3 

135.5 

1067.4 

4 

336.0 

• 

1267.9 

5 

536.5 

1468.4 

6 

737.0 

1668.9 

7 

937.4 

1869.3 

8 

1137.9 

2000.0 

9 

1338.4 

200  0 . 0 

i 

10 

1538.9 

2000.0 

M 

t it*****#**************#**#***##*******###*##*****##  ###*■***###* 

j 

THETA 

L (THETA ) 

E (THET  A ,R ) 

E(THETA,T) 

; 

0.0 

0.000 

2.3 

0.0 

151.0 

0.200 

5.6 

853.  1 

- ; 

200.5 

0.500 

5.4 

1082.9 

274.2 

0.800 

3.8 

1040.2 

<s  | 

INFINITY 

1.000 

0.0 

466.0 

***#*«*«**««***««*«**«*#*««**««««#«*#***#*««**#*«*««#«*««*«**« 


FIGURE  25 
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Appendix  I 
Glossary  of  Symbols 


a 

A 

8 

B 

X2a,v 

E0(r) 

E0(t) 

exp [x] 
h 

h0 

hl 

Hi 

k 

L(6) 

MTBF 

n 

r 

r0 

T 

0 

e 

°0 

°I 

s 

V(t) 


Producer's  risk  (Type  I error) 

Ratio  of  risks,  equal  to  (1-3)/  a 

Consumer's  risk  (Type  II  error) 

Ratio  of  risks,  equal  to  8/(l-a) 

Chi-square  variable  with  "v"  degrees  of  freedom 
with  "a"  area  above  the  variable 

Expected  number  of  failures  if  0 is  the  true  MTBF 

Expected  waiting  time  to  reach  a decision  if  0 is 
the  true  MTBF 

2.718...  raised  to  the  "x"  power 
Exponent  used  in  calculating  L ( 0 ) 

Intercept  of  accept  line  in  a sequential  test 
Intercept  of  reject  line  in  a sequential  test 
Hypothesis  concerning  alternative  i 
Ratio  of  0o  and  01#  equal  to  0o/61 

Probability  of  accepting  Hc  when  0 is  the  true  MTBF 
Mean  Time  Between  Failure 
Number  of  items  on  test 
Number  of  failures  at  time  t 

Number  of  failures  at  which  sequential  test  is  truncated 
Time  at  which  sequential  test  is  truncated 
True  MTBF 

Upper  limit  of  MTBF  (desired  MTBF) 

Lower  Limit  of  MTBF  (undesired  MTBF) 

Slope  of  accept  and  reject  lines  in  a sequential  test 
Accumulated  test  time  (in  replacement  case) 
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— X*  DISTRIBUTION  PROBABILITY  VALUES 


Values  ~ef  „ 


V 

V* 

x .995 

V* 

x .99 

V* 

x .975 

B9 

x .90 

V* 

x .80 

V* 

x .75 

V* 

x .70 

* 

1 

.0000393 

.000157 

. 000982 

.0158 

.0642 

.102 

. 148 

1 

2 

.0100 

.0201 

.0506 

.211 

.446 

.575 

.713 

2 

3 

.0717 

. 115 

.216 

.584 

1.005 

1.213 

1.424 

3 

4 

.207 

.297 

. .484 

1.064 

1.649 

1.923 

2.195 

4 

5 

.412 

.554 

.831 

1.  145 

1.610 

2.343 

2.675 

3.000 

5 

6 

.676 

.872 

1.237 

1.635 

2.204 

3.  070 

3.455 

3.828 

6 

7 

.989 

1.239 

1.690 

2.  167 

2.833 

3.822 

4.255 

4.671 

7 

8 

1.344 

1.646 

2. 180 

2.  733 

3.490 

4.594 

5.  07 1 

5.527 

8 

9 

1.735 

2.088 

2.700 

3.325 

4.168 

5.380 

5. 

6.393 

9 

10 

2.156  • 

2.558 

3.247 

3.940 

4.865 

6.179 

6.  737 

7.267 

10 

11 

2.603 

3.053 

3.816 

4.575 

5.578 

6.989 

7.584 

8.  148 

It 

12 

3.074 

3.571 

4.404 

5.226 

6.304 

7.807 

8.438 

9.034 

12 

13 

3.565 

4.  107 

5.009 

5.892 

7.042 

8.634 

9.299 

9.  926 

13 

14 

4.  075 

4.660 

5.629 

6.571 

7.790 

9.467 

10. 165 

10.821 

14 

15 

4.601 

5.229 

6.262 

7.261 

8.547 

10.307 

11.036 

11.721 

15 

16 

5. 142 

5.812 

6.908 

7.962 

9.312 

11.152 

11. 192 

12.624 

16 

17 

5.697 

6.408 

7.564 

8.672 

•10.085 

12.002 

12.792 

13.531 

17 

18 

6.265 

7.015 

8.231 

9,390 

10.865 

12.857 

13.675 

14.440 

18 

19 

6.844 

7.633 

8.907 

10.  117 

11.651 

13.716 

14.562 

15.352 

19 

20 

7.434 

8.260 

9.591 

10.851 

12.443 

14.578 

15.452 

16.266 

20 

21 

8.034 

8.897 

10.  283 

11.591 

13.240 

15.445 

16.344 

17.  182 

21  1 

22 

8.643 

9.542 

10.982 

12.  338 

14.041 

16.314 

17.240 

18. 101 

22 

23 

9.260 

10. 196 

11.688 

13.091 

14.848 

17.  187 

18. 137 

19.021 

23 

24 

9.886 

10:856 

12.401 

13.848 

15.659 

18.062 

19.  037 

19.943 

24 

25 

10.520 

11.524 

13. 120 

14.611 

16.473 

18,  940 

19.  939 

■ 20.867 

25 

26 

11. 160 

12.198 

13.844 

15.379 

17.292 

19.820 

20.843 

21,792 

26 

27 

11.808 

12.879 

14.573 

16. 151 

18,114 

20.  703 

21.  749 

22-719 

27 

28 

12.461 

13.565 

15.  308 

16.928 

18.939 

21.588 

22.657 

23.647 

28 

29 

13.  121 

14.256 

16.047 

17.708 

19.768 

22.475 

23.567 

24.577 

29 

30 

13.  787 

14.953 

16.791 

18.493 

20.599 

23.  364 

24.478 

25.508 

30 

35 

17.  156 

18.484 

20.558 

22.462 

24.812 

27.820 

29.058 

30.181 

35 

40 

20.674 

22.  142 

24.423 

26.507 

29.067 

32. 326 

33.664 

34.874 

40 

45 

24.281 

25.880 

28.356 

30.610 

33.367 

36.863 

38.294 

39.  586 

45 

50 

27.962 

29.687 

32.348 

34.762 

37.706 

41.426 

42. 944 

44.314 

50 

55 

31.708 

33.552 

36.390 

• 

38.956 

42.078 

46.011 

47.612 

49.  055 

55 

60 

35.510 

37.467 

40.474 

43.  186 

46.478 

50.614 

52.295 

53.808 

60 

65 

39.  360 

41.427 

44.595 

47.448 

50.902 

55.  233 

56.991 

58.572 

65 

70 

43.253 

45. 426 

48.750 

51.737 

55.349 

59.368 

61.698 

63.  344 

70 

75 

47. 186 

49.  460 

52.935 

56. 052 

59.815 

64.515 

66.416 

68. 125 

75 

80 

51.  153 

53.526 

57.146 

60.390 

64.299 

69.  174 

71.144 

72.913 

80 

85 

55.  151 

57.621 

61.382 

64.748 

68.799 

73.843 

75.880 

77.  707 

85 

90 

59.  179 

61.741 

65.640 

69.  124 

73.313 

78.522 

80.623 

82.508 

90 

95 

63.’  3 

65.886 

69.919 

73.518 

77.841 

83.210 

85.374 

87.314 

95 

100 

67.  312 

70.053 

74.216 

77. 928 

82.381 

87. 906 

90.  131 

92. 125 

100 

105 

71.414 

74.241 

78.530 

82.352 

86.933 

92.610 

94.894 

96.941 

105 

110 

75.536 

78. 448 

82.861 

86.  790 

91.495 

97.321 

99.663 

101.761 

110 

115 

79.679 

82.672 

87.207 

91.  240 

96.067 

102.  038 

104. 437 

106.  565 

115 

120 

83.839 

86.913 

91.567 

95.703 

100.648 

106.762 

109.216 

111.413 

120 

— X*  DISTRIBUTION  PROBABILITY  VALUES 
(continued) 

•N.  V*lue«  of  X1a^ 


V 

V* 

x .50 

V* 

* .30 

V* 

x .25 

x*.  20 

wmm 

■hHQH 

EBI 

V* 

x .025 

V* 

x .01 

V* 

x .005 

¥ 

I 

.455 

1.074 

1.323 

1.642 

2.  706 

3.841 

5.024 

6.635 

.7.879 

i 

z 

1.386 

2.408 

2.773 

3.219 

4.605 

5.991 

7.378 

9.  210 

10.597 

H 

3 

2.366 

3.665 

4. 108 

4.  642 

6.251 

7.815 

9.  348 

11.345 

12.838 

Kt 

4 

3.357 

4.678 

5.385 

5.989 

7.779 

9.488 

11. 143 

13.277 

14.860 

n 

5 

4.351 

6.  064 

6.626 

7.289 

9.  236 

11.070 

12.832 

15.086 

16. 750 

6 

5.  348 

7.231 

7.841 

8.558 

10.645 

12.592 

14.449 

16.812 

18.548 

§9 

7 

6.  346 

8.  383 

9.037 

9.803 

12.017 

14.067 

16.013 

18.475 

20.  278 

mm 

8 

7.344 

9.524 

10.219 

11.030 

13.362 

15. 507 

17.535 

20.  090 

21.955 

8 

9 

6.  343 

10.656 

11.359 

12.242 

14.684 

16.919 

19.023 

21.666 

23.589 

9 

10 

9.  342 

11.781 

12.549 

13.442 

15.987 

18.307 

20.483 

23.209 

25. 188 

11 

10.  341 

12.899 

13.701 

14.631 

17.275 

19.675 

21.920 

24.725 

26.757 

11 

12 

11.340 

14. 011 

14.845 

15.812 

18.549 

21.026 

23. 337 

26.217 

28.300 

12 

13 

12. 340 

15.  119 

15.984 

16.985 

19.812 

22.362 

24.736 

27.688 

29.819 

13 

14 

13.339 

16.222 

17. 117 

18.151 

21.064 

23.685 

26. 119 

29.  141 

31.319 

14 

15 

14.  339 

17.  322 

18.  245 

19.311 

22.307 

24.996 

27. 488 

30.578 

32.801 

15 

16 

15.338 

18.418 

19.  369 

20.  465 

23.542 

26. 296 

28.845 

• 32.000 

34.267 

16 

17 

16.338 

19.  511 

20.489 

21.615 

24.  769 

27.587 

30.  191 

33.409 

35.718 

17 

18 

17.338 

20.601 

21.605 

22.  760 

25.989 

28.869 

31.526 

34.805 

37. 156 

18 

19 

18.338 

21.689 

22.  718 

.23.900 

27.204 

30. 144 

32.852 

36.191 

38.582 

19 

20 

19.337 

22.  775 

23.828 

25.038 

28.412 

31.410 

34.170 

37.566 

39.997 

20 

21 

20.337 

23.858 

24. 935 

26. 171 

29.615 

32.671 

35.479 

38.932 

41.401 

21 

22 

21.337 

24.939 

26.039 

27.301 

30.813 

33.924 

36.781 

40.289 

42.  796 

22 

23 

22.337 

26.018 

27.  141 

28.  429 

32.007 

35. 172 

38.076 

41.638 

44. 181 

23 

24 

23.337 

27.096 

28.  241 

29.553 

33.  196 

36.415 

39.364 

42.980 

45.558 

24 

25 

24.  337 

28.172 

29.339 

30.  675 

34.382 

37.652 

40.646 

44.314 

46.928 

25 

26 

25.  336 

29. 246 

30.434 

31.795 

35.563 

38.865 

41. 923 

45.642 

48.290 

26 

27 

26.336 

30.  319 

31.528 

32.912 

36.741 

40.  113 

43. 194 

46.963 

49.645 

27 

28 

27.336 

31.391 

32,620 

34.027 

37.916 

41.337 

44.461 

48.278 

50.993 

28 

29 

28.336 

32.461 

33.711 

35.139 

39.  087 

42.557 

45.722 

49.  588 

52.336 

29 

30 

29.336 

33.530 

34.800 

36.250 

40.256 

43.773 

46.979 

50.892 

53.672 

30 

35 

34.338 

38.860 

40.221 

41.802 

46.034 

49.  798 

53.207 

57.359 

60.304 

35 

40 

39.337 

44.  166 

45.615 

47.295 

51.780 

55.755 

59. 345 

63.706 

66. 792 

45 

44.337 

49.453 

50.984 

52.757 

57.480 

61.653 

65.414 

69.971 

73. 190 

50 

49.336 

54.725 

56.333 

58.  194 

63. 141 

67.502 

71.424 

76.167 

79.512 

50 

55 

54.336 

59.  983 

61.665 

63.610 

68.770 

73.309 

77.384 

82.305 

85. 769 

55 

60 

59.336 

65.229 

66.982 

69.  006 

74.370 

79.086 

83.301 

88.391 

91.970 

* *9 

65 

64.336 

70.466 

.72.  286 

74.387 

79.946 

84.819 

89. 181 

94.433 

93. 122 

70 

69.335 

75.693 

77.578 

79.752 

85.500 

90.530 

95.027 

100.436 

75 

74.335 

80.912 

82.860 

85. 105 

91.  034 

96.216 

100.843 

106.403 

80 

79.335 

86. 124 

88. 132 

90.  446 

96.550 

101.879 

106.632 

112.338 

116.334 

; (V 

85 

84.  335 

91.  329 

93.396 

95.777 

102.  050 

107.521 

112.397 

118.244 

122.337 

85 

90 

89.335 

96. 529 

98.653 

101. 097 

107.536 

113. 145 

118. 139 

124. 125 

128.310 

90 

95 

94. 335 

101.723 

103. 902 

106.409 

113.008 

118. 751 

123.861 

129.  980 

134.257 

95 

100 

99.335 

106.911 

109. 145 

111.713 

118. 468 

124.  342 

129.  565 

135.814 

140.179 

fT’T’9: 

105 

104.  335 

112. 095 

114.381 

117.009 

123.917 

129.918 

135.250 

141.627 

146.078 

105 

110 

109.  335 

117. 275 

119.612 

112.299 

129.355 

135.480 

140.920 

147.421 

151.956 

110 

1 1 5 

114. 335 

122. 451 

124.838 

127. 581 

134.782 

141.030 

146. 574 

153.  197 

157.814 

115 

120 

119.335 

127.623 

130.059 

132.856 

140.201 

146.568 

152.215 

158. 956 

163.654 

120 

4] 


APPENDIX  III 

LISTING  OF  SEQUENTIAL 
TESTING  PROGRAM 


hhi 

h 1 f '6 . .O i OK  A ($  i f a 1 4 ) ,\r(h ) r 6 <4 ) (A ) , t H 4 ) 

1 1 

INTEGER  R,R2 

1 20 

REAL  L,L5 

13C 

DATA  A / 0.01,0.05,0.10,0.20  / 

14(1 

DATA  B / 0.01,0.05,0.10,0.20  / 

150 

PRINT  7 

160 

7 

FORMAT ( "0  ENTER  N , THETA  1 ,T0“ ) 

170 

READ ,N,O1,T0 

180 

L(1)-0.0 

190 

O(1)-0.B 

200 

PRINT  14 

210 

14 

FORMAT(“0  ENTER  0 IF  FULL  DISPLAY  IS  DESIRED** ,/,  “ ENTER 

2206 

“ 1 IF  NOT  DESIRED**) 

230 

READ  ,ICHK 

240 

IF (ICHK .NE . 0 ) GO  TO  16 

250 

II1-1 

260 

II2-4 

270 

JJ1-1 

2G0 

JJ2-4 

290 

GO  TO  17 

300 

16 

PRINT  18 

310 

18 

FORMAT  ( "0  ENTER  ALPHA, BETA  DESIRED**) 

322 

READ, AL, BE 

330 

DO  23  I» 1 ,4 

340 

IF  (AL.E3.A{I))  II1-I 

350 

IF (BE  .EQ  .B  (I  ) ) JJ1-I 

360 

23 

CONTINUE 

370 

II2-II1 

362 

JJ2-JJ1 

390 

17 

PRINT  1 

402 

1 

F0RMAT(///,24X, “SEQUENTIAL  TESTING  PLANS" ,// ,24X , 

4 12& 

“EXPONENTIAL  DISTRIBUTI0N“,//,27X,“ (WITH  REPLACEMENT)" 

42E& 

,//,24X,24 ("*")) 

432 

DO  10  I-II1.II2 

442 

DO  10  J*  J J 1 , J J2 

450 

A 1 » 1 . E-A ( I ) 

460 

BI-I.B-B(J) 

470 

AA-(1.0-B(J))/A(I) 

480 

BB«=B(J)/(1,0-A(I)) 

490 

PRINT  2,A(I),B(J),N,T0,O1 

520 

2 

F0RMAT(/,20X, “ALPHA  «“ ,F8 . 2 ,/ , 28X , “BETA  ,F8.2,/,2DX 

5106 

**N  -**,I8,/,28X,“T0  -“,F8.1,/,28X,"THETA1 

5206 

F8.1,//,2CX,16("*")) 

530 

R1«2.0*FLOAT(N)*T0/O1 

640 

DO  11  K-2,2fl0,2 

550 

X2B2RoCHIS3(B(j),K) 

560 

IF  (X282R ,LE ,R 1 ) GO  TO  12 

570 

GO  TO  13 

580 

12 

IDX-K 

590 

X21-X2B2R 

60  2 

1 1 

CONTINUE 

610 

13 

RsIDX/2 

620 

O2-2.0*FLOAT(N)*T0/CHISQ(A1,IDX) 

630 

PRINT  3,R,O0 

43 


64  Z 

3 FORMAT  ( / ,2  OX  , "R  .‘Me  , / ,2£X  ,"THETA0  «",Fe,1,//, 

C l S'  & 

2LX,16("*“)  ) 

Ci., 2 

DENOM-1 .3/01-1 .0/00 

672 

H0*-ALOG(BB)/DENOM 

,* 

66 2 

hi*alog(aa)/denom 

692 

S*AL0G(02/01)/DEN0M 

e~s  -1  <1 

/ V O 

PRINT  4,H0,H1,S 

712 

4 FORMAT (/ ,28X , "HB  -"  ,F8 . 1 , / ,28X , "H 1 «",F8.1,/f 

72  2& 

20X  , "S  *",F8.1) 

732 

PRINT  5 

740 

5 FORMAT ( / ,28X  , "PLOTTING  POINTS" ,//, 14X , "FAIL . NO.",10X, 

* 

750G 

"REJECT", 13X, "ACCEPT",/) 

760 

R2-R+1 

r 

770 

DO  15  K«1,R2 

~ 

7G0 

K1-K-1 

790 

VR  — H1+F10AT(K1)*S 

000 

VA*H0+FLOAT  ( K 1 )*S 

810 

IF(VR.LT.0.B)  VR-0.0 

023 

IF  (VA.LT.0.0)  VA-0.0 

030 

IF  (VR  . GT .T0 ) VR-TB 

04  C 

IF (VA.GT.T3)  VA-T0 

. 850 

15  PRINT  6 ,K 1 , VR ,VA 

e60 

6 FORMAT (17X,I3,11X,F8.1,8X ,F8 . 1 ) 

870 

PRINT  8 

800 

8 FORMAT (/,5X, 62 ("*" ),// ,8X , "THETA" ,8X , "L (THETA ) ",8X t 

890G 

"E(THETA,R)”,8X,"E(THETA,T)",/) 

9C  0 

0 (2 ) *0 1 

910 

0(3) “S 

920 

0 ( 4 ) *00 

930 

L ( 2 ) »B ( J ) 

• 940 

L ( 3 ) *ALOG ( AA ) / ( ALOG ( AA )-ALOG (80 ) ) 

950 

L (4 ) *A  1 

960 

DO  19  K = 1 ,4 

970 

IF  (K  »EQ  . 3)  GO  TO  21 

980 

ER (K)-(H1-L(K)*(H0+H1) )/(S-0(K)) 

993 

GO  TO  19 

1002 

21  ER(K)-H0*H1/(S*S) 

1010 

19  ET(K)-0(K)*ER(K)/FL0AT(N) 

1020 

PRINT  9,  (0(K) ,L(K) ,ER(K) ,ET(K) ,K«1,4) 

13  30 

9 FORMAT (6X ,F8 • 1 » 9X  #F5  * 3 *2 ( 1 0X  ,F8 . 1 ) ) 

1340 

L5-1.0 

105  0 

ET1*H0/FLOAT(N) 

• * 

1060 

PRINT  22 , L5  » L ( 1 ) , ET 1 

10  73 

22  FOR M AT( 6X, "INFINITY", 9X,F5. 3, 2(1 0X,F8.1)»//,5X,62( ”*"),//) 

10  00 

10  CONTINUE 

1290 

STOP 

i. 

1103 

END 

* 

1110 

FUNCTION  CHISQ (PROB  , IV ) 

1120 

P*1  . 0-PROB 

113  2 

IF ( P ) 1,4,2 

114(5 

1 PRINT  3 

■ j 

1132 

3 FORMAT ( / / , "P  IS  NOT  IN  THE  INTERVAL  (0,1),  ", 

11636 

"INCLUSIVE",///) 

| 

I; 

■ j 

/ i 
* : 

♦ i 

iV/k 

119" 
1200 
1210 
1222 
1230 
1240 
1290 
1260 
1270 
1280 
1290 
130  f*  0 
1310 
1320 
1330 
1340 
1390 
1360 


no  T n i? 

2 (6-1  . " !I  7,6,1 

4 X— B .999999E+74 

5 D=0.B 

GO  TO  12 

6 X=C.999999E+74 
GO  TO  5 

7 D-P 

IF (D-E .5)  9,9,8 

8 D-1.B-D 

9 T2«ALOG(1.0/(D*D)) 

X«T-(2!915517+0 ,802853*T +0 . 0 1 0328*T2 ) / ( 1 . 0+1 .4327BB 
*T+0 . 109269*T2+3 . 00  1308*T*T2) 

IF (P-0 .5)  10,10,11 

10  X— X 

11  X1=2.0/(9.0*FLOAT(IV)) 

CHISQ =FLOAT (IV)*(1 . 0-X 1+X*SQRT (X 1 ) ) **3 

12  RETURN 
END 
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